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ABSTRACT 

Thermal images of cold dust in the Central Molecular Zone of the Milky Way, obtained with the far-infrared 
cameras on-board the Herschel satellite, reveal a ~ 3 x 10^ M© ring of dense and cold clouds orbiting the 
Galactic Center. Using a simple toy-model, an elliptical shape having semi-major axes of 100 and 60 parsecs 
is deduced. The major axis of this 100-pc ring is inclined by about 40° with respect to the plane-of-the-sky and 
is oriented perpendicular to the major axes of the Galactic Bar. The 100-pc ring appears to trace the system of 
stable X2 orbits predicted for the barred Galactic potential. Sgr A* is displaced with respect to the geometrical 
center of symmetry of the ring. The ring is twisted and its morphology suggests a flattening-ratio of 2 for the 
Galactic potential, which is in good agreement with the bulge flattening ratio derived from the 2MASS data. 
Subject headings: Galaxy: center — ISM: clouds 



1. INTRODUCTION 

Observations of the central regions of our Milky Way 
Galaxy provide the nearest template for the study of the con- 
ditions in galactic nuclei in normal galaxies. The 3.6 x 10^ 
M© black hole at the center dominates the Central Molecular 
Zone (CMZ) which extends from Galactic longitude 1 = -1° 
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to about +1.5° and hosts the d ensest and most massive molec- 
ular clouds in the Milky Way ( [Morris & Serabyn|199 6 ). The 
CMZ may be a Galactic analog of the nuclear star-forming 
rings observed in the centers o f star forming bar red galaxies 
( Kormendy & K ennicutt 20041 I Kormendy & C omell 2004]). 
The CMZ contains^ tens of millions of solar masses oT cold 



inters tellar matter (ISM) ( Pierce-Price et aH2000,. Bally et al. 
|2010| ) and harbors several of its most acti ve sites of star for- 
mation such as Sgr A , Sgr B2, and Sgr C ('Yuse f-Zadeh et al. 
2008 Yusef-Zadeh et al. 2009 ). The gas kinematics and dis- 
tribution of near-infrared light indicate that the gravitational 
potential in the inner Galaxy is dominated by a bar with a ma- 
jor axi s inclined by about 20 to 45° with respect to our line-of 
sight ( |Binney et al. 1991 Benjamin et al. 2003 Minchev et al. 



[2007] ). 

Despite the active star formation and bar dynamics, warm 
dust h as been found t o be a relative minor constituent in the 
CMZ. |Sodroski et al.| ( p^94l ) used CQBE/D IRBE data to con- 
strain dust properties in the CMZ and Rodriguez-Fernandez 
et al. ( 2004 ) presented ISO observations. These studies found 
that 15 to 30% of the far-IR emission comes from molecu- 
lar clouds with dust temperatures of about 19 K, 70 to 75% 
arises from the HI phase with dust temperatures of 17 to 22 
K, and less than 10% is emitted by warm dust with T > 29 
K associated with the extended HII and war m (100 to 300 K) 
molecular phases traced by H3 absorption ( Oka et aL||2005 



Goto et al. 2008 ). With Herschel we can now extend tnis type 



of analysis into a new domain of spatial resolution and sensi- 
tivity. 

In this Letter, we present the first high-angular resolution 
maps of the CMZ far-IR dust emission, temperature, and col- 
umn density based o n data obtained with the Herschel satellite 
dPilbratt et al.|2010| ). 

2. OBSERVATIONS AND RESULTS 

The Galactic Center was observed with Herschel as part 
of the Hi-GA L Key-Project (the H erschel infrared Galactic 
Plane Survey, [Molinari et al. |2Q 1 Ob} . Observing strategy and 
basic data reduction are described in [Molinari et al., (2Q10a| 
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Fig. 1. — Herschel PACS 70 fim image of the Galactic Center region. Labels identify known objects that are discussed in the text. 



and Trafic ante et al. pOll ). The central 2°x2° region of the 
Galaxy was i mageTin Parallel Mod e by acquiri ng simultane- 
ously PACS ( [Poglitsch et al |201Q| ) and SPIRE ( [Griffin et al." 



Bernard et a l. ( 2010 ). Assuming constant dust properties, the 
opacity is converted into a hydrogen column density (Nh = 
Nm + Nh^) using T25o/Nh = 8.8 lO'^^ cm^/H. This value 



[2010 ) images in 5 photometric bands centered at 70, 160, 250, is estimated from latitudes b>0, and it could be higher by a 
350 and 500 jim. SPIRE was used in "bright-source mode", ^ - ^ 
to minimize non-linearities and saturation in this very bright 
region; preliminary SPIRE calibration checks indicate an ac- 
curacy within 15% compared to standard SPIRE setting (M. 
Pohlen, priv. comm.). We corrected the zero-level in our 
SPIRE and PACS images using offset values derived from the 
comparison wi th the Planck/IR AS data using the same proce- 
dure as in Bern ard et al.| ( |201()| ). 

Figure [l] shows the PACS 70 yum image. Sgr A* and Sgr 
B2 are the brightest spots, and the entire area is filled with 
intense far-IR emission punctuated by several InfraRed Dark 
Clouds (IRDCs) seen in silhouette. A large bubble dominates 
at 70 Jim from 0° ^ / ^ -h0.20° but disappears at longer wave- 
lengths (see SPIRE 250 yum in fig. At 250 yum Sgr B2 is the 
brightest feature (it is saturated at z50 yum). Two bent chains 
of emission extend above and below the mid-plane from Sgr 
B2 to Sgr A*, where they appear to cross, and continue on 
to Sgr C. This structure can a lso be identified in the SCUBA 
map ( Pierce-P rice et al.|2 000 ), but only the unique multi-band 
far-IR/submm capabilities of the Herschel PACS and SPIRE 
cameras make possible the determination of its temperature 
and column density. 

Opacity and dust temperature are estimated from a pixel- 
to-pixel fit to the 70-350 jim data (the resolution is matched 

to the 350 fim images 25'') using the DustEM model 

( [Compiegne et al.|[2011 ) following the method described in 



factor 2 - 4 in the Galactic plane (Ber nard et al.| | 2010| ) and 
could change our mass estimates accordingly. The tempera- 
ture and column density maps (Figs. |3]and|4]) represent the 
integrated contribution of dust along tne line of sight. High 
opacity, cold dust appears in silhouette against the warmer 
background even at 70 yum (fig. [ij. Although the derived tem- 
perature might be underestimated, a temperature map without 
the 70 yum data provides similar results. 

The detailed study of the great bubble visible in the 70 yum 
and the temperature maps (figs. |T]andj3]) will be presented in a 
forthcoming paper. This Letter is dedicated to the analysis of 
the large-scale filaments visible in the far-IR (fig|2]), that trace 
a continuous chain of cold and dense clumps (Td^st ^20 K; 
N(H) > 2 X 10^^ cm"^) organized along an cxi-shaped feature 
that dominates the image between the Sgr C complex at / = 
359.4° and Sgr B2 at 0.7°. The total projected extent is -1.4°, 
or about 180 pc, f or a solar Galactocentric distance of 8.4 kpc 
IReid et al.|2009 l. Warmer dust with relatively lower column 
ensity (TjL? >25 K, N(H) < 2 x 10^^ cm'^) fills the interior 
of the oo-shaped feature. 

3. A 100-PARSEC ELLIPTICAL RING OF COLD AND DENSE 
MOLECULAR CLOUDS: THE X2 ORBITS 

The oo-shaped feature in the temperature map in fig. |3]c( 
incides with high-column density material (N(H) > 2 x 10' 
cm"^) organized in a continuous chain of irregular clumps 
(see fig. |4]). Most of the cold, high column density clumps 
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Fig. 2. — Herschel SPIRE 250 jim image of the Galactic Center region. 



in the [/", b~] and [/^, h^] arms of this feature are seen in sil- 
houette against the 70 fim background (fig. [T]); the same is not 
true for clumps along the [/^, and [/", d"] arms. We con- 
clude that the [/", b~] and [/^, b^] arms are in front of the bulk 
of the warmer dust emission (i.e. between us and the Galactic 
Center), while the [/", b~] and [/^, b^] arms are located in the 
background. 

We use the CS survey by |Tsuboi erar| ( |1999| ) (their fig. 4) 
to extract the radial velocity of CS condensations positionally 
coincident with the cold and dense peaks on varius positions 
along the oo-shaped feature in figs. [3] and |4| The position- 
velocity information is averaged in latitudebins of 1.5', so 
the information can be extracted with some degree of approx- 
imation but should be adequate to represent the mean velocity 
of the various clouds and clumps. The full span of velocities 
found toward each position is reported in fig. |4] 

The [/", b^] arm shows steadily increasing (m absolute val- 
ues) negative velocities from the center of symmetry toward 
Sgr C indicating approaching material with increasing rela- 
tive speed. Following the [/",/?"] arm from Sgr C back to 
the centre we find negative velocities decreasing and chang- 
ing to positive (receding) values at / -359.8°. Continuing 
along the arm the radial velocities are still positive 

and slowly increasing, indicating material receding at higher 
and higher speeds up to Sgr B2. Finally along the [l^ ,b~] 
arm we see velocities decrease again going back toward the 
center of symmetry of the cx)-shaped feature. The distribution 
of radial velocities along the oo-shaped feature seems to indi- 
cate an ordered, rotating pattern around its projected center of 
symmetry. 



3.1. A simple model for the oo-shaped feature: a 100-pc 
elliptical twisted ring 

A simplified toy model of the oo-shaped feature is built as- 
suming that the material is distributed along an elliptical orbit 
(projected onto the Galactic Plane) with semi-axes [a, b] and 
major-axis position angle Op. For simplicity, a constant orbital 
speed Norb is assumed. An additional sinusoidal vertical oscil- 
lation component with a vertical frequency and a phase 6^ 
is added to the equations of the ellipse describing the orbit. 
Adopting a coordinate system with the x axis oriented toward 
negative longitudes, y axis pointing away along the line of 
sight from the Sun, and z axes toward positive latitudes (see 
Fig. [5]), the position and radial velocity of each point along the 
orbit IS described as a function of the polar angle 6t (counter- 
clockwise from the x axes) as 



( X = a cos 6t cos 6p - b sin 6t sin 6^ 
\ y = a cos 6t sin Op -\- b sin 6t cos 6^ 
' z = Zo sinv^iOp - 0^) 
. = -Vorb cosiOp -h Ot) 



We let a, b, Op and Vorb vary until they match the observed 
projected morphology of the ellipse and the distribution of 
radial velocities at 20 positions along the oo-shaped feature. 
We keep zo fixed at 15 pc, which is half of the measured b- 
extent of the oo-shaped feature; we also fix y^=2 in units of 
the o rbita l frequency, as indicated by the projected twist shape 
(see § 33 ). The fitting procedure minimizes a pseudo-;^^ value 
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(Vr - ypbs? 
^^obs 



(1) 
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Fig. 3. — Temperature map of the Galactic Center region. The Log-color stretch extends from 15K to 40K on Sgr A*. 



where Av^z,^ is half of the velocity range toward each of the 
20 test positions (see fig. |4]). 

The best fit is obtainea for an ellipse with a = lOOpc, 
b =60pc, Op = -40°, and Vt,rZ?=80 kms"^ 6^ ~0 implying 
that the ring crosses the midplane along the X axes. These 
numbers should be taken with caution, as the velocities are 
derived by eye from a plot. However, the discrepancy between 
the model radial velocity and the observed ones is within ±25 
km s"^ The only exception is the 50 km s~ | clo ud (see fig. (sj 
located close in projection to Sgr A* (see |3.2| ). 

The CMZ has been the subject of many studies in mil- 
limeter molecular tracers, and although the interpretation is 
difficult for the extreme confusi on along the line of sight. 



several fea tures were identifie d. Sofue (1995]) analysed the 
CO data in |Bally et aL ( |1987|) and iden t ified t wo large-scale 
structures (the Arms' ). |Tsuboi et al.| (|1999 ) recognized a 
large-scale feature in their CS position- velocity diagram, that 
they call "Galactic Center Bow" which matches the [/",/?"]- 
[/^, b^] arms of the cx3-shaped feature. The cold, high column- 
density dust features identified in the Herschel images, used 
as a proxy for the densest molecular gas, now indicate that the 
previously identified "Arms" and the "Galactic Center Bow" 
are parts of a single, twisted and elliptical ring rotating around 
Sgr A*. Figure [5] presents a sketch of the proposed arrange- 
ment of what we will call thelOO-pc ring. 

A rough estimate of the ring mass is obtained by integrat- 
ing the N(H) map within a contour surrounding the 100-pc 
ring which encloses material with N(H) 2 x 10^^ cm"^ (for 
a corresponding solid angle of ~ 4 x 10"^ ^r). The contri- 
bution of background or foreground gas is removed by sub- 
tracting a column density of N(H) = 2 x 10^-^ cm~^, the mean 



N(H) value in the interior of oo-shaped feature. The resulting 
ring mass is ~ 3 x 10^ M© (it is ~ 4 x 10'^ M© without back- 
ground subtraction). Within the uncertainties, the derived ring 
mass agrees with 5.3x10^ M© estimated from SCUBA maps 
dPierce-Price et al.|2000| ) and 3x10^ M© from CO data ( |Dah-1 
|menetal.|1998| ). 

The 100-pc ring traces gas and dust moving about the nu- 
cleus on an elliptical orbit whose major axis (the straight red 
line in the right panel of fig. [5]) is perpendicular to the major 
axis of the great Galactic bar, oriented 20 to 45° with respect 
to our line-of-sight (see fig. |5]). The 100-pc ring is likely lo- 
cated on the so-called orbit system, which precesses at 
the same rate as the solid-body angular velocity of the G alac- 



tic Bar, r esulting in stable, non-intersecting trajectories ( Bin 
ney et al.|19 91 ). These X2 orbits are enclosed within another 
larger system of stable orbits, called xi, elongated along the 
Galactic Bar. The prominent Sgr B2 and Sgr C star-forming 
complexes are located at the two projected longitude extrema 
of the 100-pc ring where gas on xi and X2 orbits may collide, 
producing strong shocks which may trigger the formation of 
massive molecular cloud formation as well as enhanced star 
formation r ates by the sho ck focusing mechanism suggested 



by Kenney & Lord| ( [1991j) for the spiral-bar interface of M83. 

It is tempting to speculate that the 100-pc ring constitutes 
the remnant of a more persistent dusty torus that may have 
played a role in past AGN phases of our Galaxy. The recent 
Fermi-LAT detection of a large gamma- ray bubble emanating 
from the Galact ic Center ( Su et al. 2010,ICrocker|201 Ul Zubo^ 
vas et al.||2011| ) provides possible evidence for such past ac- 



tivity. Mid-IR interferometry of AGNs reveals co mpact dusty 
structures with a radius of a few parsecs, (e.g. |Jaff'e et aL 
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Fig. 4. — Atomic hydrogen column density map of the Galactic Center region. Color scale shown is logarithmic and extends from 4 x 10^^ cm ^ in the d ark 
regions, to 4 x cm 2 in the brightest region corresponding to Sgr B2. Velocity information is extracted from CS spectroscopic cubes jTsuboi etal.|l999t for 
the gas counterparts positionally associated with the dense dust clumps; the range values indicate the range where emission is observed. 



|2QQ4[ IRadomski et al.|2QQ8] ), but there is additional evidence 
from Spitzer imaging and spectroscopy suggesting that larger 
and more persistent str uctures play a ro le in shaping o bserve d 
AGN properties (e.g. |Shi e t al.""2006| [Quillen et al. 2006). 
Model predictions (e.g. IKrolik 2007 ) for torus column den- 
sity in excess of 10^"^ cm"^, and orbital speed -100 kms"\ 
seem to be consistent with the parameters we estimate for the 
100-pc ring. 

3.2. The relation between the 100-pc ring and Sgr A*" 

Sgr A*is noticeably displaced -10' (24 pc in projection) 
towards the negative longitude side of the geometrical center 
of symmetry of the 100-pc ring. The ring displacement with 
respect to the supposed location of the Galactic Center is con- 
sistent with the overall displacement of dense gas within the 
CMZ. Two-thirds to three-quarters of the cold interstellar gas 
and dust within the inner few hundred parsec s of the nucleus is 
at pos itive longitudes with respect to Sgr A*(|Bally et al.|1988| 
Morris & Serabyn 1996). A similar asymmetry is present in 
the observed radial velocity on the 100-pc ring close to the 
location of Sgr B2, as opposed to Sgr C (see fig. |4]). This may 
indicate an additional pattern motion of - 10-20 kms"^ of 
the 100-pc ring as a whole along its major axes and toward 
us, implying a 'sloshing' motion back and forth that could be 
explained as the respo nse of the gas to an m=l mode in the 
Galactic bar potential ( [Morris & Sera byn 1996 ). 

The projected proximity of the 20 and 50 kms~^ clouds to 
Sgr A*, has led to the suggestion (e.g. Hermstein & Ho|20"Q5 



of the ring than to its back; Sgr A* may then not be in the ring 
center, nor lie along the ellipses axes. Proximity to Sgr A* 
could be the source of the roughly 40 km s"^ deviation of the 
radial velocity of the 50 km s"^ cloud from the value predicted 
by the model of the rotating elliptical ring. 

Given the extent of the 100-pc ring, its orbital motion is 
driven by the enclosed gas and bulge stars. If, however, 
ring material passes in proximity of Sgr A* the gravitational 
pull due to the strong mass concentration can be substan- 
tial. Indeed, the 3.6 x 10^ M© black hole in Sgr A* domi- 
nates the gravitational potential only to a radius of a few par- 
secs. Beyond that distance the potential is dominated by the 
stellar mass in the bulge with a radial density profile of the 
form p(r) ^ Po(r/riy)~^ whe re a ^ rp ~ 0.34pc and 

po ~ 2.1 X 10^ Mo/pc^ (Trip pe et aLl|2008| ). We then treat 
the gravitational pull from the Sgr A"*" mass concentration as 
an additional component over the general rotational motion. 
Using Newton's 2"*^ law: 



GMc(r) 



Av 

^0 



(2) 



that they are physically interacting with Sgr A*. If they are 
part of the 100-pc ring, then Sgr A* must be closer to the front 



where r is the distance between Sgr A* and the two molec- 
ular clouds (that is assumed similar for the two clouds if they 
belong to the 100-pc ring), Mc(r) is the mass enclosed in a 
sphere of radius r centered on Sgr A* and to is the timescale 
for the cloud passing in front of Sgr A*. The latter can be 
written as to = Ax/Vorbr where Ax is the projected distance 
between the 50 and 20 kms"^ clouds and Vorb is their orbital 
speed. Eq. |2]can then be rewritten 
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Fig. 5. — Sketch of the proposed 3D structure and placement of the 100-pc ring. Left panel: the ring is represented by the thick color line (red and green mark 
positions above and below ^7=0°. The line of sight to the Sun is along the Y axes. The thin full and dashed lines represent the major and minor axes of the ellipse. 
Right panel: Top view of the 100-pc ring with the proposed location of major clouds. The thick black dashed lines represent the innermost x\ orbits. The position 
of SgrA* along the line of sight is the one corresponding to the distance it should have from the front portion of the 100-pc ring to justify the velocity difference 
between the 20 and 50 km s~^ clouds if due to the gravitational pull of the mass concentration around Sgr A* (see text). 



\GMc{r)^x 

\ AV Worb 



(3) 



Using this formulation to estimate Mc(r), the measured 
projected distance between the nearest edges of the two 
clouds of Ax ~ 7pc, and assuminsan orbital speed of 
80km s"^ from our toy model fit, Eq.^ provides r -20 pc. 
In other words, to explain the diff'erence in radial velocities 
between the 20 and 50 km s"^ clouds would require that Sgr 
A* is closer than 20pc to the foregreound arm of the 100-pc 
ring. 

3.3. The ring's twist: constraints on the flattening of the 
Galactic potential 

Perhaps the most striking aspect of the 100-pc ring is its 
twist. Apparently the gas in the ring oscillates twice about 
the mid-plane each time it orbits the nucleus once. Assum- 
ing the elliptical axes and the orbital speed derived from the 
toy model, the orbital period is ~ 3 x 10^ years; the pe- 
riod of the vertical oscillations is therefore P^ ~ 1.5 x 10^ 
years. The coherence of the ring's morphology suggests sta- 
bility over a few revolutions. Thus, the relationship between 
Vo and P^ can be used to infer constraints on the flattening of 
the gravitational potential in the CMZ. Assuming as a gross 
approximation that the matter (gas and stars) is distributed in 
a uniform-density po slab of thickness h, P^ at a distance R 
can then be expressed as 



^Gpoh 



(4) 



A perturbation on a particle rotating in the flattened slab 
potential will result, in a P^-order approximation (Binney &] 



|Tremaine|1987| ), in a harmonic oscillation both in radial and 
vertical directions. The vertical oscillation frequency can be 
expressed as 



(5) 



R,0 



The gravitational field along the z direction of the slab (ap- 
proximated as infinite) is 



Fz = InpoGh 
so that the period of the vertical oscillations is 



In 



(6) 



(7) 



As noted above, the oo shape that the ring displays in pro- 
jection implies that 



or, using Eqs. |4]and[7j 



h. 



R 

2 



(8) 



(9) 



The shape of the 100-pc ring would therefore imply a flat- 
tened potential in the CMZ that is consistent with a x/z = 0.5 
axial ratio bulge mass distribution, remarkably similar to the 



value of 0.5-0.6 recently derived ( [Rodriguez-Fernandez & 
Combes| [M08 ) fitting the 2MASS data to a potential model 
that includes the bulge and the nuclear bar. 
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4. CONCLUSIONS 

The far-infrared multi-band images provided by the Hi- 
GAL survey using the Herschel PACS and SPIRE cameras 
provide evidence for a 100x60 parsecs elliptical ring of cold 
and dense clouds orbiting around the Galactic Center on X2 
orbits elongated along the minor axis of the Galactic cen- 
ter Bar. The ring is twisted with a vertical frequency that is 
twice the orbit frequency, indicating a flattening-ratio of 2 for 
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